Noisy kink in microtubules 
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ABSTRACT. We study the power spectrum of a class of noise effects generated by means 
of a digital-like disorder in the traveling variable of the conjectured Ginzburg-Landau-Montroll 
kink excitations moving along the walls of the microtubules. We have found a 1/P noise with 
a G (1.82-2.04) on the time scales we have considered. 
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One can find extensive, descriptive presentations of microtubules (MTs) in many biological 
papers. Here, we shall give only elementary definitions as follows. They are ubiquitous protein 
polymers in eukaryotic cells belonging to the category of biological filaments and making the 
most part of the cytoskeleton. They are hollow cylinders 25 nm in outer diameter and 17 nm 
in the inner one, with lengths ranging from nm to mm in some neuronal cells. The walls of 
the cylinders are usually made of 13 (the seventh Fibonacci number) protofilaments laterally 
associated. The surface structure of MT walls is very interesting Structurally, MTs are quasi 
one-dimensional chains of tubulin polar dimers (negative a and positive (5 monomers, each of 
4 nm in length) undergoing conformational changes induced by the guanosine triphosphate 
to diphosphate (GTP-GDP) hydrolysis. The whole assembly process of a MT is due to the 
hydrolyzation of GTP. The cation Mg ++ is essential to increase the affinity of tubulin for binding 
GTP and thus for generating MTs. Moreover, a unique dynamical property is the so-called 
dynamical instability Q which is a random growth and shrinkage of the more active plus ends 
of MTs. It has raised much interest in recent years |J. 

An energy-transfer mechanism in MTs by means of Ginzburg-Landau-Montroll (GLM) kin- 
klike protofilament excitations has been discussed by Sataric, Tuszyhski and Zakula § in 1993. 
Also sine-Gordon (SG) solitons have been discussed by Chou, Zhang and Maggiora in 1994 ||. 
We recall that various types of solitons have found interesting applications in biological physics 
(Davydov's model ||, DNA/RNA |J7||). Usually, in order to get nonlinear differential equations 
one performs a continuum limit for some lattice models in which discreteness effects are ne- 
glected. Some authors have shown that such effects might be important and suggested various 
ways of including them in continuous differential equations ||. The digital disorder we shall 
comment on and use next is just a possible way to incorporate discreteness in a moving solitonic 
pattern. 

In 1993, Rosu and Canessa || introduced a digital-like disorder in the Davydov /3-kink 
leading to the 1/P noise, with a ~ 1, in the dynamics of that kink, and also commented on 
the multifractal features of the dynamics of the /3-kink. The procedure is as follows. In the 
traveling variable £ = x — vxt of the slowly moving kink (the acoustic Lorentz factor ~ 1) 
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one puts Xi = xq + Axj, where xq is the position of the center of mass of the kink, or its central 
position along the chain, and Ax are small random displacements around xq (say, in the interval 
±1). In the calculations, one can fix io = 0. For each random position Axi chosen from a 
uniform distribution, one calculates by means of a fast-Fourier-transform algorithm the noise 
power spectrum of the time series of the signal (considered to be the kink) in order to get the 
time correlations of the fluctuations of the signal, i.e., 

;(l [ K(0e 2mft dt\ 2 ^ , (1) 

where 0<t<r~l// and K(£) is the kink function, and the brackets stand for averaging 
over ensembles. This approach to noise effects has been taken from the literature on the self- 



organized criticality (SOC) paradigm, see, e.g., |10[. On the other hand, the standard treatment 



of noise effects when they originate in thermal fluctuations is by means of the Langevin equation 
method. Vails and Lust [|ll]] have studied the effect of thermal noise on the front propagation in 
the GL case. They have found a crossover between constant-velocity propagation at early times 
and diffusive behavior at late times. 

The purpose of the present work is to investigate the noise produced by the same type of 
disorder as in |J in the case of the GLM kink conjectured in MTs. 

The main assumption in M is that the assembly of tubulin dimers/dipoles (T> n ) form a quasi 
one-dimensional ferrodistortive system for which the double-well on-site potential model is a 
standard framework 

V(p n ) = - l ~AVl+ l l BVi. (2) 

The variable T> has been identified in H with the amount of /3-state distorsion vertically projected 
(the /3-state is defined as having the mobile electron within the /3-monomer). Moreover, in Q 
a GL hamiltonian/free-energy with intrinsic electric field and dissipation effects included led to 
the dimer Euler-Lagrange dimensionless equation of motion (EOM) in the traveling coordinate 
of the anharmonic oscillator form (with linear friction) 

\A\V" - -yav K v' - F(V) = , (3) 

where 7 is the friction coefficient, a = \A\ r )\lMv 2 soun( ]. and F(V) = AV - BT> 3 + qE, with q 



denoting the effective charge of a single dimer of mass M, and E the magnitude of the intrinsic 
electric field. This EOM is known to have a unique kink solution given by the formula [12] 

K({) = " + tt^S) s ■ + 7=2 {1 ~ ta " h(/5f/2)) • <4> 

where K = T> / y/\A\/ B is a rescaled dipole variable, (3 = (b — a)/y/2, whereas a and b are two of 
the solutions of the cubic equation 

F{K) = {K- a){K - b){K - c) = K 3 - K - a = , (5) 



where a = q ^rj^- E and no units are used, where uq = y/\A\/B ~ 1.4- 1CP 11 m is the amplitude 
of the dimer displacement (shift of the double- well potential). Notice that the GLM kink is thin. 
Its width is wk = -s ~ 0.7no- 

As we said, the type of digital disorder we consider here is very close to the ideas of the SOC 
paradigm that we understand in the broad sense of both spatial and temporal scaling of the 
dynamical state of the system [^]. The spatial scaling (self-similarity) is of the (multi) fractal 
type while the temporal scaling leads to l/f a noises. 

A strong motivation for dealing with digital dynamics is the possibility of generating broken 
symmetries and therefore of having various types of dynamical phase transitions Thus, 
digital disorder, though might look a rather ad-hoc approach, focuses on both self-organized 
properties of MTs, i.e., to driven steady states with long-range spatio-temporal correlations, 
and to (dynamical) phase transitions, since digital dynamics allows for symmetry breaking. 

Our results are displayed in Figures 1 and 2 and show that the noise introduced by the digital 
disorder in the GLM kink variable is practically of the 1 /f 2 (Brownian) type on the time scales 
we have considered. At low frequencies there is the known cross-over to a white noise due to the 
finite system size, which is moving to lower frequencies as the size of the system is increased | jlq| . 
On the other hand, the deviation from the power law at high frequencies is an artifact due to 
the so-called aliasing 16]. The Brownian noise we have obtained is not unexpected since it is a 



common occurence in SOC models in any finite dimension [17], and only mean- field calculations 



reproduced the 1// noise [18 |. 
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Finally, we would like to mention that if polarization does not exactly follow the displace- 
ment one needs a system of two coupled partial differential equations leading to two coupled 
traveling kink waves. Of course, the nonlinear models seem to be too simple-minded for the 
MT complexity. Nevertheless they provide guides for further insight and perhaps some partial 
answers. 
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Fig. 1 

Double logarithmic plot of the power spectrum of the digital noise perturbing the motion of the GLM 
kink moving along MTs at constant n k = 2m/s. The fitted slopes are as follows: (a) -2.032 (b) -2.0175 
(c) -2.0164 for the time scales corresponding to that of the dimer, ten times bigger, and hundred times 
bigger, respectively. The errors in the slopes are at the level of 0.0003 for each case. 
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Fig. 2 

The same plot as in Fig. 1 for v K = lOOm/s. The fitted slopes are (a) -2.0083 (b) -1 .8641 (c) -1.8205 
(d) -2.0000 for temporal scales of 10 -2 , 10 _1 , 10°, and 10 1 times that of the tubulin dimer, respectively. 

The level of the errors is the same as in Fig. 1. 
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